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C$ ' Abstract 

We address the important issue of stabilizing the dilaton in the context of superstring cosmology. Scalar 
potentials which arise out of gaugino condensates in string models are generally exponential in nature. In a 
cosmological setting this allows for the existence of quasi scaling solutions, in which the energy density of the 
scalar field can, for a period, become a fixed fraction of the background density, due to the friction of the 
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background expansion. Eventually the field can be trapped in the minimum of its potential as it leaves the 



scaling regime. We investigate this possibility in various gaugino condensation models and show that stable 
ON • solutions for the dilaton are far more common than one would have naively thought. 
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I. INTRODUCTION 



Scalar fields in cosmology have been extensively studied over the past few years. One of the most intriguing 
areas in which they occur is Superstring theory, where the presence of the dilaton field is vital. Its vacuum 
expectation value (VEV) determines both the gauge and gravitational coupling constants of the low energy 
theory and also fixes the scale of supersymmetry (SUSY) breaking through the gravitino mass, m^n. Therefore 
realistic models require a VEV of order one (in Planck units), and TO3/2 ~ 1 TeV. 

Unfortunately, in string theory the dilaton potential is flat to all orders in perturbation theory, which of course 
means there is no way of obtaining a stable VEV for the field. This problem has to be overcome through some 
nonperturbative effect. The most promising possibility is through the formation of condensates of gaugino fields 
at an energy scale of around 10 14 GeV [Q. The resulting scalar potential for the dilaton is then a combination of 
exponentials and polynomials in the field. A detailed investigation of these condensate models has demonstrated 
the need for at least two condensates to form if the dilaton potential is to develop a minimum at a realistic 
value although with a negative cosmological constant (these are the so-called "racetrack" models) . 

An alternative proposal has recently been suggested as a method of obtaining a minimum for the dilaton field, 
and it has the advantage of relying on only one gaugino condensate ||,|| . In this scenario the Kahler potential 
(which determines the kinetic terms of the dilaton in the action) requires string inspired nonperturbative cor- 
rections. A detailed analysis of these models || indicates that it is possible to have a minimum with zero or 
small positive cosmological constant. One additional positive feature that emerges is that the nonperturbative 
corrections can lead to a solution of the "moduli problem" for the dilaton || (fields with masses in the TeV 
range but which decay so slowly that they spoil nucleosynthesis), by giving it a huge mass. 

Although attractive, both kind of models still have several problems associated with them. One is the difficulty 
of achieving inflation which was carefully studied by Brustein and Steinhardt Q a few years ago. Taking a 
model of multiple gaugino condensation (as those studied in Q ) they argued that the kinetic energy associated 
with the dilaton field would dominate over its potential energy until <fi (the canonically normalized field, related 
to the usual dilaton by ReS* = e^) would settle near a minimum of its potential. This obviously excludes 
inflation from happening, at least with the dilaton as the inflaton field, leaving the possibility of <f> settling down 
to a minimum and then inflation being driven by other fields. However this second possibility also presented 
serious problems as the models they studied had a negative cosmological constant. Nonperturbative corrections 
to the Kahler potential can cure this latter problem, but in both cases the potentials are exponentially steep 
in the strong coupling regime. This would lead us to expect the dilaton to roll past the minimum rather than 
acquiring its VEV, which seems to be a major problem that superstring cosmology needs to address. 

In this paper we turn our attention to the possibility that other matter fields rather than the dilaton drive 
the evolution of the Universe. Recent attention in cosmology has turned to the investigation of scaling solutions 
in models with exponential scalar field potentials [p|-|lO[ . These models are of particular interest because if the 
background dynamics are dominated by some matter source other than the field itself (i.e. radiation, dust, 
vacuum energy) then it is possible for the field to enter a scaling regime as it evolves down its potential. In this 
regime the friction term from the expansion of the Universe balances the kinetic energy of the field allowing it 
to enter this scaling era. Attractor solutions exist HJl(| where the energy density in the field becomes a fixed 
fraction of the total energy density. 
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This intriguing behaviour can be applied to the case of the dilaton field arising from string theory. Under the 
assumption that it is evolving from somewhere in the strong coupling regime of its potential we show how, for 
a wide range of initial field values in the presence of a background dominated by a barotropic fluid, the dilaton 
enters a quasi scaling regime as it evolves down the potential, inspite of its steepness. This scaling behaviour 
eventually ends as the field enters the minimum of its potential, by which time it has slowed down sufficiently 
for it to simply oscillate about it, losing energy and eventually becoming fixed with a realistic VEV. 

In section II we introduce the concept of scaling solutions with exponential potentials. In section III we 
demonstrate how this can be successfully adapted to both the racetrack and modified Kahler potential models 
of gaugino condensates. Solutions are presented analytically and numerically showing how the dilaton field is 
stabilized in its minimum. We conclude in section IV. 



II. SCALING SOLUTIONS WITH EXPONENTIAL POTENTIALS 

In this section we will describe some general features concerning the cosmological evolution of scalar fields with 
exponential potentials. Let us consider a scalar field <f> with a potential energy density given by V = Vb e ~ AK ^, 
with k 2 = 8ttG and A and Vo constants, which is evolving in a Friedmann-Robertson- Walker (FRW) Universe 
containing a fluid with barotropic equation of state p 7 = (7 — l)p 7 , where 7 is a constant (0 < 7 < 2, for 
instance 7 = 4/3 for radiation domination or 7 = 1 for matter domination). The equations of motion for a 
spatially flat FRW model with Hubble parameter H are 

k 2 

H = - — {p 1 + Pl + 4> 2 ) 

p 7 = -3H(p 7 + p 7 ) (I) 

• dV 
4> = -3H(f>-— , 



subject to the constraint 



H 2 = ^{p 1 + \<P 2 + V) . (2) 



In the previous equations we have assumed that the only interaction between <j> and the other matter fields is 
gravitational. In what follows we set k 2 = 1 but it can easily be reinstated. For exponential potentials the 
asymptotic behaviour of this system can be obtained analytically, and we shall focus on the solution for very 
steep potentials, namely A 2 > 37, as these are of most interest in string theory. Let us first of all review the 
structure of the solutions which, as is well known ||, contain a late- time attractor solution. For that purpose 
it is useful to proceed as in |nj, defining the variables x = 4>/\/&H and y = W /V3H and using the logarithm 
of the scale factor TV = ln(a) as the time variable. The previous system Eq. (Q) becomes 

x' = -3x + A^|y 2 + \x[2x 2 + 7 (1 ~ x 2 ~ y 2 )] 

V =-X\[\xy+ Z -y[2x 2 +1 {l-x 2 -y 2 )] (3) 

H' = -^H[2x 2 +1 (l-x 2 ~y 2 )} , 

where a prime denotes a derivative with respect to N. In terms of these variables, the constraint Eq. @ becomes 
x 2 + y 2 + p 7 /3H 2 — 1 or, in other words, for p 7 > we have the bounds < x 2 + y 2 < 1. 
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FIG. 1. 

Plot of the evolution of (f> versus N = ln(a). The scalar potential is given by V = e~ A ^ with A = 85, 7 = 1 and the initial 
conditions are </>o = 0.5, <j>o = 0.2 and Hq = 1. 



The evolution of the held in terms of N can be seen in Fig. 1, for A = 85, Vq — 1 and initial conditions 
(pa = 0.5, cf>o = 0.2, and Hq = 1. After an initial increase driven by the initital velocity of the held, we see how 
the friction term (or, in other words, the expansion of the Universe) dominates and freezes (f> at a constant value 
for a considerable amount of time, until the held reaches a scaling regime corresponding to the critical points 
x c = ^/3/27/A and y c = ^3(2 — "f)"f/2X 2 This behaviour can also be obtained by solving the system of 



equations Eq. (jjj) under the following assumptions: 



• Stage I (pre scaling regime): Since A > 1, and H is decreasing, y is very small in this early stage of the 
evolution and can be neglected in Eq. (|^). The equation for x now becomes 



-3a; + -x[2x 2 +7(1 ~ x 2 )] , 



and its solution is 



x = 1 



2 ^0^3(2-7)^ 

n fc- 



-1/2 



(4) 



(5) 



where xq is the initial condition for x (at N = 0). The solution for <j) can be obtained by integrating 
Eq. (||) (recall that x — cj>' /y/E), and is given by 



3(2- 7 ) 



sinh 



— sinh 



-3(2- 7 )AT/2 



(6) 



where <Pq is the initial value of the held. As N increases this solution tends to a constant value <po given 
by 
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mfi±^l . (7) 



Obviously for zero initial velocity, 4>o will reduce to 

• Stage II (scaling regime): as mentioned above, the scaling regime is defined by constant (i.e. critical) 
values for x and y. We then obtain H — i7 e -37Ar / 2 and 



where now we have an explicit dependence on the characteristics of the potential (i.e., A and Vo). 

It is important to note that the background evolution is being determined by the additional matter fields 
present, given by 7. It is not the <fr field which is driving the evolution. The fact that the simultaneous evolution 
of the scalar field together with the background causes the former to reach a scaling regime, inspite of the 
steepness of its (exponential) potential, suggests that this could be also the case for the kind of potentials 
arising from SUSY breaking via gaugino condensation, where the superpotential depends exponentially on the 
dilaton field. In the next section we will apply these solutions to the gaugino condensation models studied in 
I and §. 

III. SCALING SOLUTIONS WITH GAUGINO CONDENSATES 

A. Two condensate potentials 

Multiple gaugino condensation (or racetrack) models have been extensively studied in the literature 11 3] . 
Essentially the idea is to consider a strong type interaction in the hidden sector of our effective supergravity 
(SUGRA) theory, which is governed by a nonsemisimple gauge group. The superpotential of such models will 
then be expressed in terms of a sum of exponentials which conspire to generate a local minimum for the dilaton. 
To be more precise, the scalar potential in any N=l SUGRA model |l3] is given by 



V = e K \W\' 



Mr 1 + ^ 

W I V J w 



(9) 



where K is the Kahler potential, W is the superpotential and the subindices i, j represent derivatives of these 
two functions with respect to the different fields. Given that we are interested in superstring derived models, 
and in particular in studying the hidden sector of the theory, both K and W will be dependent on the dilaton 
(S) and the moduli (Tj, i = 1,2,3) fields. In fact we know that in the case of orbifold compactifications the 
tree-level Kahler potential is given by: 

3 

K = -\og(S + S)-Y / ^g(T l +f l ) , (10) 

j=i 

and we will restrict our study to the case of a hidden sector interaction governed by two gauge groups, SU(iVi) 
x SU(-/\?2) under which we have Mi[N\ +N±) and M2(N2 + N2) "quark" representations with Yukawa couplings 
to a set of singlet fields. For simplicity, we will assume an overall modulus T = Ti = T 2 = T 3 and a generic 
singlet field for each of the gauge groups, A\ and A2 respectively. In this case, the superpotential is given by 
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,M 4 / % _ a(S 



■Mi 



(ii) 



where a; = 8?^/^, ft = 2(3A^ - Mi)/N h d t = iV i (327r 2 e)( M< / JV *- :L ), jj(T) is the Dedekind function and the self 
coupling of the Ai field is set to hi = 1. [Note iV, is not the same as TV, the number of e- foldings defined earlier]. 

These kind of models were thoroughly studied in j|] , so let us summarize their main features: the presence of 
the r\ function, imposed by the requirement of target space modular symmetry and, in general, the T-dependence 
of the potential ensures the presence of a minimum for T ~ 1.2 (in Planck units) ]l3| ], independently of the 
particular gauge groups and/or matter representations (provided the dilaton acquires a VEV); also there exist 
minima in the ImS" direction if both condensates have opposite phases. Finally it was shown the existence 
of many examples for which there is a minimum in the ReS* direction at the phenomenologically acceptable 
value ReS" ~ 2 (remember that ReS = gj t ri„ g ), with a reasonable (~ 1 TeV) gravitino mass but always with a 
negative value of the potential energy. A typical example is shown in Fig. 2(a), which perfectly illustrates the 
problem that Brustein and Steinhardt pointed out in their paper: the steepness of the dilaton potential, which 
would prevent the field from settling down at its (negative) minimum, instead allowing it to run over the tiny 
maximum towards infinity. 
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FIG. 2. 

(a) Solid line: plot of the scalar potential V (in logarithmic units) vs ReS 1 for two condensates with gauge groups SU(6) X 
SU(7) with Mi = 2 and M2 = 8 matter representations repectively; dashed line: our exponential approximation given by Eq. (|l2|). 
(b) Evolution of Re S vs N for the two condensate potential plotted in (a) with Ho = 1, 7 = 4/3, the different initial positions 
being ReSo = 0.6 (dot-dashes), ReSo = 1 (dots), ReSo = 1.65 (dashes), ReSo = 2 (solid) and the initial velocities given by 
ReSo = ReS'o/4. The thick solid line represents the solution for the scaling regime Eq. (|l8|), for the exponential potential of 
Eq. @. 



We would like to study the evolution of the dilaton field (more precisely of its real part) before it settles to the 
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minimum. We have already seen that a single exponential scalar potential can be solved analytically leading to 
scaling solutions. However there are major differences between this and the more realistic gaugino condensate 
models. First, the dilaton couples not only gravitationally but also directly to the matter fields. For simplicity, 
and to avoid making any assumptions about specific models, we will neglect this effect throughout this paper. 
Second, in the case of two condensates, the superpotential Eq. contains two different exponentials of ReS 1 , 
therefore the scalar potential Eq. (^) will have all the different terms coming from \ W\ 2 . In particular there will 
be a mixture of exponential and polynomial terms in the dilaton. And finally, the differential equation for <f> 
written in Eq. ([!]) is meant to be obeyed by a canonically normalised field, and ReS* is not. As it is well known, 
in the SUGRA Lagrangian the kinetic terms for scalar fields are given by K^D^iD^^ , which in the case of 
KeS introduces an extra factor of l/(2RcS') 2 . Therefore, the correct procedure in order to study the evolution 
of the dilaton in an expanding Universe would be to solve the system in terms of the canonically normalised 
field 0, with ReS* = or, alternatively, to modify Eq. ([!]) accordingly in order to account for the non canonical 
kinetic terms. 

Fortunately these problems can be overcome. First we note that, even though the two condensates in Eq. ( |i"l| ) 
have to be carefully fine-tuned to produce a minimum for Re5 at the right value, for most of the evolution 
towards such a minimum only one of them (the one with a smaller on value) will dominate. For simplicity, we 
also keep the matter fields constant at their minimum value Ai m [ n = (3jVj/ 'diMi)~ 2 ' exp(— 2ai5 , m ; n / /3j) /rj(T) 2 
during the evolution of HeS. Therefore in the region we are studying, the superpotential can be approximated 
by, for example, the first condensate; moreover it is also easy to show that, in Eq. (^|), the term proportional 
to dW/dAi dominates among those within the brackets. In conclusion, our scalar potential can be very well 
represented in the region before the minimum by the following expression 

2 



V = 



min 



\/6^i(2ReT)7 ? (r)' 3 i 



which corresponds to considering only the dominant term mentioned above and setting ReS 1 = 1 everywhere 
but in the exponential. The result of such an approximation can be seen in Fig. 2(a), represented by the dashed 
line (with, as everywhere in this article, ReT = 1.2), and it is good enough to justify our use of Eq. ( |l2| ) when 
studying the evolution of Re5* away from the minimum. 

Concerning the second problem, that of Re5 not being a canonically normalised field, we start by showing 
the exact numerical result of the evolution of the <f> field, plotted in Fig. 2(b) in terms of ReS = versus N, for 
a radiation dominated Universe (i.e. 7 = 4/3) and initial conditions Ho — 1 and ReSo = Rc5o/4. The different 
lines correspond to different initial conditions for Re5. It is remarkable how the behaviour in the first stages of 
the evolution is very similar to that shown in Fig. 1 for a pure exponential potential, the difference appearing 
after N ~ 11 e- foldings and due to the presence of a minimum in this case. Depending on the initial position 
of the field it may or may not fall in the minimum, but what is clear is that if the field reaches the scaling 
regime the former will certainly happen, and that occurs, as we can see, for a very wide range of values of ReSo, 
contradicting the general belief that only for a very narrow range of initial values, all around the minimum, 
would the dilaton settle at its minimum. In fact the top curve of Fig. 2(b), which is the only one that does not 
end up at the minimum, corresponds to an asymptotic value of the field Re<So = beyond the maximum of 
the potential. 
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We have checked that the scaling solution shown in Fig. 2(b) is very well represented by Eq. (g) in the case of 
the potential being given by Eq. (Jl|), i.e. A = 2a t and V Q = {dxMj^A^l'^^ 1 /^/6Ni(2ReT)rj(T)^\ 2 . That 
is, as if ReS were a canonically normalised field. This approximate solution corresponds to the thick solid line 
shown in Fig. 2(b). 

We turn our attention to explaining this scaling behaviour of the non canonical dilaton field by trying to solve 
Eq. ([!]) analytically for this case. In terms of the normalised dilaton, </>, the equations to solve are 

x' lfi = -3x + Xe^^yj + ^x [2x^ + 7(1 - x\ - y\)\ 

y'4 = -Ae^yfaw + \vA^\ + 7(1 " x l ~ Vl)\ (!3) 
H' = ~H[2xl + 1 {l-xl-yl)] , 

where x§ and y^ are the x and y of Eq. (^|). Note that the presence of a more involved potential requires the 
replacement A — > Ae*. We can proceed as in the case of the pure exponential and solve for the two different 
stages defined before. Solving for Stage I is trivial, as the key point in this regime is to neglect any dependence 
on the potential and therefore the solution (for ip) is identical to that in Eq. (^|). Then 

/ / 2 \ -^ /(2 ^ } 

ReS I = Re5 ^ e~i^ N + Jl + -^ e -^-^ , (14) 

V V^4 V 1 ~< J 

where ReSo = e^°, with <^o given by Eq. (^). However solving for the scaling regime of Stage II is much 
more complicated. To begin with, the form of the potential and the range of values of <f> we are interested in 
guarantees that x,y 1, and therefore we can solve for H, obtaining the usual result H = Hoe~ 3lN ^ 2 , once 
again indicating that the background fields are determining the evolution of the Universe. As for the other two 
equations, let us rewrite them in terms of ReS 1 and a new variable xs = ReS" /V6 (note that ys = 2/0) ■ Again 
in the same approximation of small x, y we have for the two first equations 

3 ,„ s . , / 3 2 2 



2 (2- 7 )x s +Ay-(ReS)'i/S 

/3 3 
^x s y s + --yy s ■ (15) 

Obviously, for the scaling regime observed in Fig. 2(b) (i.e. Stage II) we must have x' s — y' s = 0. The second 
equation in Eq. (|l5| ) will give us then the expected solution for xs, analogous to the pure exponential case 

%l ■ 

whereas from the first equation we find that 

2 _ 3(2~ 7 ) 7 1 
Vs ~ 2A2 (ReS) 2 ' { ' 

that is, ys does not seem to reach a critical value but instead has a dependence on (ReS*) -2 . However, given 
the size of A (> 20) and the range of values we are considering for ReS* (between 0.3 and 2), the deviation of ys 
from its expected critical value, given by y c s = 3(2 — 7)7/(2A 2 ), is not going to be significant, as it is obvious 
from Fig. 2(b). In any case let us compute this correction which modifies the pure scaling result by a factor of 
e(N). Our ansatz is then 
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ReS » = T N+ J ln [ 9Hl(2- lh ) + ^ ■ (18) 

Substituting into Eq. ( |l7| ) and using the definition of ys (= \/Voe ~$ RcS / (y/3H)) and the solution for H we 
obtain 



e(N) =~ln 
A 



3 Win' 2V °* 



(19) 



which is indeed a very small numerical correction to the standard result (in fact x' s = e"(N) ~ 0), and the only 
noticeable effect of having solved for a scalar field with an exponential potential but non minimal kinetic terms, 
a very encouraging result. 

The pure exponential approximation will eventually break down as the dilaton approaches its minimum, 
R-eSmm. As can be seen from Fig. 2(b), if the field is in its scaling regime, it will not have enough energy to 
go over the maximum of the potential. Instead, it will oscillate around the minimum with an exponentially 
damped amplitude, settling down quickly to its final value. A simple estimation of the number of e-foldings 
needed to reach the minimum can then be obtained from the scaling solution Eq. ( |l8| ) by equating ReSn to 
R-eS'min. One obtains, 



iVmin = ^~ 

37 



AReSVnin — In 



2V n A 



2 



(20) 



,9tf 2 (2- 7 ) 7/ 

(where we have ignored the e correction as this result is accurate enough). Therefore once we have defined the 
example we are working with and the background, we will have N m i n . In fact, with this result we can also 
calculate H mm 



which, remarkably enough, is independent of Hq. 

We have estimated N m - ln and H mm in a radiation dominated Universe (7 = 4/3) for a number of hidden sector 
gauge groups with HeS ~ 2 and m 3 / 2 ~ 1 TcV, obtaining the almost invariant result 

N min ~ 11 

H min ~ 5.1CT 10 M P . (22) 

Assuming a radiation dominated Universe, where H oc T 2 /Mp, this implies T m in ~ 10 13 GeV. The invariance 
of this result can be more or less understood by rewriting Eq. ( p0| ) in terms of the gravitino mass and the gauge 
group factors of the two condensates. The requirement of the model to be phcnomcnologically viable fixes most 
of those parameters leading to a constant value for Af m ; n and, therefore, H m i n and T m j n . 

From these results it is possible to estimate the range of initial conditions that will eventually lead to a stable 
dilaton. As long as the field enters the scaling solution before reaching its minimum, we can ensure that it will be 
stabilized. For this to happen, a suficient condition is to take an initial value Re5o between the scaling solution 
at N — (i.e. the constant term in Eq. (|l8|)) and ReS^m, and an initial velocity such that the asymptotic 
solution KeSo in Eq. (|l4|), is smaller than ReS'min. Namely we obtain the following bounds for ReSo and Xq, 

1 / 2U A 2 \ 
A ln ( v 9go 2 (2°-7)7 j <Re5o<Regmin 
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(23) 

Some examples can be seen in Fig. 2(b). Note that this region gets smaller for decreasing Ho, as expected. 
There will also be a few initial conditions outside these limits that will still lead to a stable dilaton, such as 
initial values between the maximum and the minimum, values close to the lower bound on RoSo or allowing 
for negative initial velocities. Still, even if we do not take these into account, it is clear that there is a sizable 
region in parameter space that allows the dilaton to evolve to its minimum and stay there. 

A further possibility is to consider an inflationary scenario. This will correspond to the case of a very small 
(and changing) 7. Note that it is the barotropic fluid that is producing inflation, so that we are not considering 
any kind of dilaton driven inflation. We can see from Eq. ( |T^ ) that for 7 ~ the scaling solution is practically 
horizontal. It would take a large amount of e-foldings for this solution to reach the minimum of the potential. 
Furthermore, for some models of inflation, 7 would be rapidly changing and one does not expect the field to 
follow the scaling solution exactly in these cases pp| . Nevertheless, a few e-foldings of inflation can open up a 
large region of parameter space. Unless the energy density is completely dominated by the dilaton, the almost 
constant H of an inflationary scenario will quickly freeze the field at a constant value that will then lead to 



a scaling solution during reheating (7 = 1). The bound for Xq in Eq. (23) still applies, but is maximized for 
7 = 0. One then expects most of the region of parameter space to evolve to a stabilized dilaton if there is an 
initial small period of inflation. 



B. One condensate with Knp 

In this section we will perform the same analysis of the cosmological evolution of the dilaton field for a 
different class of models proposed more recently || ||] . These consist of a single condensate which gets stabilized 
by the presence of nonperturbative corrections to the Kahler potential, which form has been suggested in ]l4| ]. 
Therefore we are dealing now with a scalar potential given again by Eq. (^), where the superpotential is simply 

W = -!^e~ aRcS ', (24) 

with a = 8ir 2 /N for a SU(N) group and C — — N/(32ir 2 e). The Kahler potential is now more involved, 
K = Ko + K np , where Ko is defined by Eq. ( |Io|) (with an overall modulus) and the nonperturbative correction 
is parameterised as 

^ np = _^=log(l + e- s (^-^)) , (25) 

that is in terms of three constants So, D, and B, the first of which just determines the value of ReS 1 at the 
minimum. Therefore this description is effectively made in terms of only D and B, which are positive numbers. 
It has been shown || that for a wide range of values of both parameters it is possible to generate a minimum 
at Sq with zero cosmological constant. The shape of these potentials is again very similar to that of the two 
condensate models. 

Another feature of this ansazt for the nonperturbative corrections is that it is very well approximated in the 
region ReS 1 < So by the following expression: 
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and therefore the scalar potential is given by 



V 



-D(l-y/So/B*S) 

2ReS 



4(1 - D^/So/RcS + 2aReS) 2 



\W\ 



(27) 



4 + 3Dy/S /ReS 

In this second example, even though we have a single condensate to start with, the dependence of the potential 
upon ReS is much more complicated, as can be expected from the form of the nonperturbative corrections to 
K, given by Eq. (p6|). To be able to write Eq. (|27]) in the usual exponential form V = Vbe _ARcS we set ReS" = 1 
everywhere except in the exponents for which we use a linear fit. We obtain 



V 



,2(Dy/So~ + 2a) 2 ( C 



-(2a+D/So)RcS 



(28) 



(4 + 3ZV3o) \v(T) 6 . 

that is a pure exponential, where the exponent depends on the value of D. As it was mentioned before, for a 
given hidden gauge group, there exist a series of values of (£>, B) for which V has a minimum at Sq with zero 
cosmological constant. That means that the cosmological evolution of a particular hidden sector interaction will 
be different depending on which values of the pair (D, B) we are considering. This can be clearly appreciated 
in Fig. 3 where we plot the cosmological evolution of the dilaton field versus the number of e-foldings N for the 
same gauge group, SU(5), and initial conditions Ho = 1, 7 = 1, ReSo = 0.15, ReSo = 0.6. 
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FIG. 3. 

Evolution of Re S vs N for the potential of Eq. ( g7j ) with gauge group SU(5) and initial conditions Ho = 1, 7 = 1, Re5o = 0.15, 
ReSo = 0.6. The different lines correspond to different values for D in Eq. (p^): D = 1 (dash-dotted), D = 3 (dotted), D = 10 
(solid) and D = 20 (dashed). 



Depending on the value of D (the corresponding B is fixed in order to have a zero cosmological constant at 
the minimum given by Sq) wc sec that the field will fall into the minimum provided that this minimum is not 
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too fine-tuned. This is determined by the magnitude of D, as was described in j5j, the smaller it is the bigger the 
value of B (and the amount of fine-tuning in the potential) is to set V — at the minimum. Therefore in this 
particular case we see how for D < 1 the minimum becomes too fine-tuned and is unable to stop the field from 
rolling past the maximum. Unfortunately these small D solutions are precisely the ones which correspond to a 
largest hierarchy between the dilaton and the gravitino masses ||, and therefore provide us with a solution to 
the "moduli problem" for this field. However, for the range of values of D which give a satisfactory stabilization 
of the dilaton we can generate ratios between those two masses of up to 300 which would, in most cases, be more 
than enough to solve the above-mentioned problem. Changing the value of 7 would correspond to a change in 
the minimal value of D from which onwards the dilaton would fall in the minimum. That is, a bigger (smaller) 
value of 7 corresponds to a bigger (smaller) value of the minimal D, and therefore to a smaller (bigger) hierarchy 
between the dilaton and gravitino masses. 

A similar analysis to the one performed in the previous section for the two condensates model can be done here 



concerning the analytical solutions of the evolution equations (for which we use Eq. (28)). Solving for Stage I 
is identical as before (see Eq. (|l4|)), as the system of differential equations given by Eq. (|J) does not depend on 
the characteristics of the potential in this regime. For Stage II we can obtain the analogous of Eq. ( [l8"| ) and in 
this case it is defined by A = 2a + D/S and V = e D 2(D v >% + 2a) 2 C 2 /[(A + iD^/S^-q 12 ]. Expressions for the 
number of e- foldings needed to reach the minimum, A m j n as well as for H m { ni are also obtained by replacing 
in Eqs. (|2^J^) the current expressions for a and Vq. Once again, the fact that we are imposing a consistent 
phenomenology (i.e., ReS* ~ 2 and m 3 / 2 ~ 1 TeV) at low energies implies that the numerical values of N m i n 
and i? m in are very similar to those obtained in the previous section. In particular we can express A m ; n as 



-2 

A min = — 
37 



\Mp J \ ^(4 + 3L>V5£) J V 2 ^ 



(29) 



For the example shown in Fig. 3, namely 7=1, H — 1, we obtain, to a very good approximation, the almost 
invariant value A m j n ~ 19 which can be translated into H m \ n ~ 10~ 13 . 



IV. CONCLUSIONS 

We have examined the issue of the cosmological evolution of the dilaton field in gaugino condensation models 
of SUSY breaking. By studying the behaviour of a scalar field <f> with an exponential potential, we have been able 
to obtain analytic expressions for its time evolution in the presence of the Hubble parameter and a dominating 
background consisting of other matter fields. It turns out that under such circumstances the field (f> tends to 
enter a scaling regime defined by constant values of both the potential and the velocity of the field relative to 
the expansion of the Universe, provided that this potential is steep enough. 

Encouraged by these promising results we turned to apply them to the dilaton field in two particular models 
of SUSY breaking: multiple gaugino condensation (or racetrack) and one condensate with nonperturbative 
corrections to the Kahler potential. However two major differences arise with respect to the ideal case we had 
dealt with before: the shape of these potentials is not exactly exponential and the dilaton field is not canonically 
normalised. We have shown that these two problems can be easily overcome, and we have found very accurate 
analytic approximations to explain our numerical results. In both cases it is possible, contrary to the general 
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belief, to stabilize the dilaton at its minimum for a large range of initial conditions (i.e., initial values for its 
position and velocity with 7 ~ 1) despite of the steepness of the potentials in their strong coupling regime. A 
previous small period of inflation (i.e., 7 ~ 0) opens up even more the region of parameter space leading to a 
stable dilaton. Moreover, the number of e-foldings required to do so and the value of the Hubble parameter at 
the minimum seem to be fixed by the requirement of having a successful phenomenology at low energies. 
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